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QED on a momentum lattice * 

J.B. Kogut and J.-F. Lagae ^ 

Department of Physics, University of Illinois, 1110 West Green Street, Urbana IL 61801, USA 

We investigate the possibility of doing momentum space lattice simulations as an alternative to the conventional 
method. The procedure is introduced and tested for quenched QED2 and quenched QED3. Interesting physical 
applications to unquenched QED3 and quenched QED4 are also briefly discussed. 
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1. Introduction 

The basic idea of the momentum space lattice 
method is very straigthforward: one just stud- 
ies continuum QED on the torus (or more gener- 
ally on any compact manifold) . Fourier expansion 
then leads to a discretized problem and the num- 
ber of degrees of freedom can be made finite by 
imposing an ultraviolet cut-off. For simplicity, we 
will consider here that the length L of the torus 
(size of the box) is the same in each direction. 
However, finite temperature type situations can 
also easily be handled. 

After gauge fixing (to the Feynman gauge) and 
Fourier expansion, the Euclidian action of QED 
in d dimensions can be recast in the form: 
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where, instead of the continuum fields 7p and A^j_, 
we have used the following dimensionless vari- 
ables: 



.(g) - eL^^-^^A^iq) 



and (3 — L^'^^^'>/e'^ p = niL. In the above for- 
mula, the q^' s take integer values and the fc^'s 
are integers or half-integers depending on whether 
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the boundary conditions for the fermions are pe- 
riodic or antiperiodic. On an N"^ lattice (= defini- 
tion) the sums over k and k' are truncated in such 
a way that N positive frequency modes are in- 
cluded in each direction, q then takes all the pos- 
sible values of the momentum transfer between k 
and k' . 

The method has two main advantages: 1) the 
coupling of the fermions to the gauge fields is non- 
compact and 2) there is no fermion doubling. The 
first is crucial for physical applications to QED3 
and QED4 as will be discussed below. The second 
is simply a consequence of the fact that what we 
consider is nothing more than a particular formu- 
lation of continuum QED. The method also has 
two potential problems, which we should discuss 
here. They are: 1) The fact that the fermionic 
matrix is no longer sparse in momentun space 
(due to the convolution with the gauge field) 
and 2) The breaking of gauge invariance induced 
by imposing an ultraviolet cut-off. The first is 
easily dealt with by using repeated Fast Fourier 
Transform in matrix multiplications. This brings 
the computer time requirements from 0{N^) to 
0{N In N), which is still ~ IniV slower than the 
position space method; however, because of the 
non-compact character of the coupling mentioned 
above, smaller lattices are generally sufficient in 
momentum space, so that on the overall, we can 
still gain. The problem of gauge invariance is well 
known from conventional perturbation theory in 
quantum electrodynamics: 11^;^ is not transverse 
when computed with an ultraviolet cut-off. The 
cure to this problem is to use a gauge invariant 
regularization: for example, Pauli-Villars. In the 
case of QED2 and QED3, the procedure simply 
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consists in dividing the usual fermionic determi- 
nant by the determinant of a heavy fermion of 
mass M and taking the hmit A — > oo (number of 
modes going to infinity) before letting M — > cx). 
In 4 dimensions, the procedure is slightly more 
involved and would include two regularizing de- 
terminants jl). Coming back to the previous case, 
the implementation of the algorithm is straight- 
forward: in the Hybrid Monte Carlo algorithm, 
one replaces the usual pseudo-fermion potential 

by (j)* ■^Mi-^m^\n)~^ ^M't'- ^^^^ ^ computcr 
time is concerned, the additional multiplication 
by Am doesn't make much difference. 



2. QED2 
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from this approximation remains interesting. We 
measured both (V"/") and the mass gap on lattices 
of various "physical size" (eL=20,30,40) for vari- 
ous number of modes (2 x 2, 4x4, 6x6, lOx 10) and 
for various values of the bare fermion mass (rang- 
ing from m/e — 0.02 to m/e — 0.2). Some of our 
results for {tp^j) are compared on fig. 1 with those 
obtained from the position space non-compact 
formulation of QED2. Although in both cases the 
extrapolation to an infinite ultraviolet cut-off ap- 
pears rather difficult, an agreement between the 
two methods seems possible. It would certainly 
be interesting to investigate this in greater detail 
particularly in relation with the claim that the 
chiral condensate should diverge in the quenched 
Schwinger model (remember that we have only 
the zero topological sector here). We have also 
measured the mass gap in this model. The results 
are presented in Fig. 2 in the form of M{p)/\p\ 
where M is the dynamical fermion mass and p 
is the lowest momentum available on the lattice 
(po ~ ii/L,pi = 0). As one should expect, the 
results converge faster with increasing number of 
modes on the smaller lattices. However, the sig- 
nal for chiral symmetry breaking becomes better 
and better as the physical size of the lattice is 
increased. Extremely good signals for chiral sym- 
metry breaking, are in fact a general feature of 
the momentum space lattice method, as will be 
verified below. 



Figure 1. A comparison between position and 
momentum space measurements of the chiral con- 
densate. eL = 20 and m/e ranges from 0.02 to 
0.2 (from top to bottom) 



Because it is computationally cheap, QED2 has 
been considered as a test case. We restricted our 
attention to the zero topological sector of the the- 
ory, i.e. configurations for which the fiux of the 
electric field through the torus is zero (whereas 
in general 4> = 27m). This restriction would 
be a drastic step in the full Schwinger model 
(where the chiral condensate comes from non- 
trivial topology 1^). However, in the quenched 
case to be discussed here, the theory resulting 



3. QED3 

As is well known, the breaking of chiral sym- 
metry in unquenched or finite temperature QED3 
relies entirely on weak long range interactions. 
As a consequence, conventional lattice studies of 
this problem have been plagued by strong finite 
size effects and quantitative measurements of 
the critical temperature or a determination of the 
phase diagram as a function of the number of fla- 
vors have remained beyond reach. The difficulty 
essentially comes from the fact that since large 
lattices are needed and the number of sites is lim- 
ited by computer requirements, the lattice spac- 
ing a is necessarily large and one falls inevitably 
into the strong coupling regime (remember that 
/3~^ = e^a in (2-|-l) dimensions). The interest of 
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Figure 2. Dynamical versus bare fermion mass 
on 2 X 2 and 6x6 lattices of size eL=40(top), 
30(middle) and 20(bottom) 



applying the momentum space lattice method to 
this problem, comes from the fact that, in this 
formulation, the coupling of the fermions to the 
gauge field is non-compact. Therefore, one can 
increase the "physical size" of the lattice while 
remaining relatively close to the continuum limit. 
Our preliminary results so far mainly concern the 
quenched case where we have checked that the 
method gives a very clear signal for chiral sym- 
metry breaking : In the relevant range, the dy- 
namical mass at zero momentum is always much 
larger than the bare fermion mass. 

4. QED4 

One of the most striking features of the numer- 
ical simulations of quenched non-compact QED4 
is that the results do not agree with the analytical 
prediction of Miransky. It has been argued that 
this is a consequence of induced 4-fermi interac- 
tions on the lattice and shown that the critical ex- 
ponent are expected to vary continously along the 
critical line in the (a, G) plane However, until 
now, no procedure has been found for modifying 
the lattice action in such a way that the results 
of the simulation would move (along the critical 
line) in the direction of the Miransky-Bardeen- 
Leung-Love point. We would like to argue here 
that the momentum space lattice method is a 



good candidate for performing this task. The 
first step consists in identifying the origin of the 
4-fermi interactions on conventional lattices. An 
obvious candidate is revealed by performing ex- 
actly the integration over the gauge field in the 
path integral (which is straightforward for the 
non-compact action after gauge fixing) . The pro- 
cedure gives rise to an effective action which con- 
tains a 4-fermi interaction, a "tadpole improved" 
kinetic energy term, the usual current-current 
interaction and higher order interactions. The 
derivation shows clearly that the 4-fermi term is 
a direct consequence of the compact character of 
the coupling between fermions and gauge fields 
and would therefore not be there on a momentum 
lattice. As a consequence, we expect that the mo- 
mentum space method will give results which are 
much closer to those of Miransky or at least will 
allow us to check whether these results survive 
the introduction of vertex corrections. 

5. Conclusion 

We have seen that the procedure of momentum 
space lattice simulation is practical and gives a 
very clear signal for chiral symmetry breaking for 
quenched QED both in 2 and 3 dimensions. We 
hope to be able to report soon on physical ap- 
plications to quenched QED4. The unquenched 
theory in lower dimensions is also currently being 
tested. 
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